Current-driven and field-driven domain walls at nonzero temperature 
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We present a model for the dynamics of current- and field-driven domain-wall lines at nonzero 
temperature. We compute thermally-averaged drift velocities from the Fokker-Planck equation that 
describes the nonzero-temperature dynamics of the domain wall. As special limits of this general 
description, we describe rigid domain walls as well as vortex domain walls. In these limits, we 
determine also depinning times of the domain wall from an extrinsic pinning potential. We compare 
our theory with previous theoretical and experimental work. 
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I. INTRODUCTION 

Current-driven domain wall motion was first predicted 
and observed by Berger in the eighties fi* 2 " It was not until 
the discovery, in the nineties, of the spin transfer torque 
mechanism, 3 '- that research on current-driven domain 
walls took off. Spin transfer torques on a domain wall 
can be understood on an intuitive level: the electrons 
which constitute the current have spin, and this spin ro- 
tates when it passes through the domain wall, as it aligns 
with the domain wall magnetization. By conservation of 
spin, there is an opposite torque on the magnetization 
of the domain wall, which leads to a net displacement 
of the domain wall in the same direction as the elec- 
tric current. Later, to explain some discrepancies with 
experiments, a so-called dissipative spin transfer torque 
(sometimes referred to as the non-adiabatic spin transfer 
torque) was added to the model^^ The value of the di- 
mensionless parameter /?, which gives the strength of this 
torque, has been the subject of much debate. By now it 
is generally accepted that (3 is of the same order as a, the 
Gilbert damping parameter, but not necessarily equal to 
it A———— Furthermore, neither f3 nor a needs to be con- 
stant. They depend on the properties of the material and 
are most likely also temperature dependent. 

Several properties of current-driven domain walls have 
been studied. One particular subject of interest is the 
velocity of the domain wall. The effect of an exter- 
nal pinning potential and the depinning behavior, both 
without and with thermal fluctuations, were investigated 
exp er iment allyi 1 2 ' 1 3 1 14 There is also interest in more com- 
plex, higher dimensional domain wall models, like vor- 
tex walfsj 15 ' 16 i 17 These are especially attractive from an 
experimental point of view since their dynamics, such 
as precession in a potential and transformations be- 
tween vortex walls and transverse walls, can be directly 
observed^^^ 

With a few exceptions ; 21 ' 22 ! 23 ! 24 ? 25 most theoretical pa- 
pers are restricted to the zero-temperature case. This is 
unfortunate, since experiments on current-driven domain 
walls are usually done at room temperature, and the rela- 
tively large current will heat up the sample even further. 
Furthermore, especially in the presence of a pinning po- 



tential, thermal fluctuations are anything but negligible. 
Domain pinning can be necessary, for instance to pre- 
cisely locate a domain wall, but thermal depinning can 
also be useful to lower the critical current. This means 
that it is very important to understand the influence of 
thermal fluctuations on the behavior of the domain wall 
precisely and thoroughly. Here, we present a unified pic- 
ture of previous work on domain wall motion at nonzero 
temperature that involved two of us, as well as new re- 
sults following from this unified picture. 

We start this paper with the Landau-Lifschitz-Gilbert 
equation, including both the reactive and the dissipative 
spin transfer torques. In Section [TTJ we apply a varia- 
tional principle to derive the general equations of motion 
for a domain-wall line including thermal fluctuations. In 
Section IIIII we investigate the velocity of current-driven 
and field-driven rigid domain walls at nonzero tempera- 
ture, while in Section HVl we look at vortex domain walls 
in more detail. Both the rigid and the vortex domain 
wall are special cases of the general description given 
in Section [TTJ In both sections, we start with the zero- 
temperature case, and then investigate the influence of 
thermal fluctuations. Pinning potentials are included in 
both models, and we examine thermal depinning. Each 
section is divided into three subsections: one in which the 
model is described, one in which we present our results, 
and a final one to make a comparison with other work. 



II. DOMAIN- WALL LINES 

In this section, we derive the equations of motion for a 
domain- wall line in a more detailed way than in previous 
work 21 and present new results obtained from this model. 
We derive the Fokker-Planck equation of the system, and 
determine the stochastic behavior under the influence of 
temperature. The model is then considered in the ab- 
sence of extrinsic pinning so that there is only intrinsic 
pinning due to magnetic anisotropy. We determine the 
drift velocity of the domain wall as a function of the cur- 
rent through the system. Finally, we compare our model 
with other theoretical and experimental work available 
in the literature. 
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A. Model 

Magnetization dynamics including spin-transfer 
torques are described by the Landau-Lifschitz-Gilbert 
equation^ 



dt 



w s -v n = n x h 
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where CI is a unit vector in the direction of the magne- 
tization. In this expression, the first term on the r.h.s. 
contains a contribution due to the effective field, which 
is written as the functional derivative of the micromag- 
netic energy functional of the system and an external 
magnetic field H(x, t) — —SEMM[^(x,t)]/HSSl(x,t) + 
gH ex t(x,t)/H. For clarity, we denote positions in three- 
dimensional space with an arrow, and the directions of 
the magnetization by bold symbols. The second term 
describes Gilbert damping, which is characterized by the 
dimensionless parameter a. The term proportional to v s 
on the l.h.s. is the reactive spin-transfer torque. The 
term proportional to (3v s on the r.h.s. corresponds to the 
dissipative spin-transfer torque and is characterized by 
the dimensionless parameter (3. The velocity v s is given 
by v s = —a PJ c /\e\, where a is the lattice constant, J c is 
the charge current and P is its spin polarization. To fa- 
cilitate a variational approach, we note that the Landau- 
Lifschitz-Gilbert equation is obtained from 
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where R is a dissipation functional, and where S = So + 
Sdrivc denotes the action of the system. The dot is a time 
derivative ft = dd/dt. The action for the magnetization 
dynamics in the absence of current and field is given by 

5 [n(af J t)]= fdtf ^JftA(n(^)).^M 



m(x, t) ■ v 2 n(f , t) - Kj_n 2 y {x, t) + k z q 2 z (x, t) 



dt 



HA(n(x, t)) 



dSljx, t) 
dt 



E M mMx, t)} 
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In this expression, J is the spin stiffness, K± > and 
K z > are the hard- and easy-axis anisotropy, respec- 
tively. The function A(fi) is the vector potential of a 
magnetic monopole which obeys Vfj x A(O) = and is 
required to reproduce the precessional motion of fi(af, t) 
around the effective field. 

The external field B oxt (x, t) and the reactive 
spin-transfer torque are determined from the action 
Sdrive[^(^)*)]- We take the external magnetic field 
to be spatially homogeneous and time independent 



B ex t(:zf, t) = B cxt . The action is given by 

Sdrivc[fi(£, *)] = J dt J^-(gB ext -Sl(x,t) 

+A(0)- (*v V)fl(x,*)j , (4) 

with g a positive constant. The dissipation functional 
that describes the dissipative spin-transfer torque and 
the Gilbert damping is written as 



dt 



v s -v n{x,t) 



(5) 



As the last ingredient, we take into account thermal 
fluctuations. We add to the effective field in Eq. ([1]) 
stochastic contributions h such that H — > H + h, where 
h has white-noise correlations 

(hi(x,t)hj(x ,t')) =ai-j5(x - x)5(t - t') , (6) 

(h i (x,t))=0, (7) 

and where the indices i,j € {x,y,z} label Cartesian co- 
ordinates. The strength cr.y is given by the fluctuation- 
dissipation theorem cry = 5ij2ak^,Ta? /h, which assures 
that, in the absence of field and current, the Boltzmann 
equilibrium distribution 



-E MM [n]/k B T 



(8) 



is reached after sufficiently long times. In principle, 
Eqs. (J2HS1) describe the full magnetization dynamics. Ob- 
taining results on finite-temperature average drift veloc- 
ities of domain walls is however very cumbersome, espe- 
cially in the presence of extrinsic pinning. We therefore 
use a variational method. 

We obtain the specific form of our variational ansatz 
by varying the action in Eq. (J3]) for a time-independent 
magnetization. Using i7 = (sin 9 cos 4>, sin 9 sin </>, cos 9) 
we find 9" = (K z /J) sin 2 9^ where the primes de- 
note derivatives with respect to x. This equation has 
domain- wall solutions tan(0/2) = exp{±[a; — A^]/A} and 
4> G {0,7r}, where A = \J J / K z is the domain- wall width 
and X is the position of the domain wall. 

The variational ansatzwe use is tan(#dw/2) = exp{[x— 
X(z,t)]/X} and <f> = 4>${z,t). A domain wall is now de- 
scribed by two collective coordinates X(z,t) and 4>o(z,t) 
that represent the position of the domain wall line and 
the chirality at this position, respectively. Note that z 
is the coordinate along the line, and also that there are 
other possibilities for the exact form of the ansatz, such 
as choosing a different domain- wall charge. 

We choose the magnetic field pointing in the positive z 
direction B ext = B z z, B z > 0. Furthermore, the current 
is taken in the positive x direction. We find that the 
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action in terms of the collective coordinates X(z,t) and 
4>o{z, t) is then given by 



S{X,fo] = 



K 



J {X 
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where N = 2XL y L z /a 3 is the number of magnetic mo- 
ments in a domain wall, with L y ,L z the length of the 
sample in the y and z direction, respectively. We now 
substitute the ansatz in Eq. ((5]) to obtain the dissipation 
functional as a function of the collective coordinates 



R[X, <j> ] = N 



ah 



dt 



dz 




(10) 



Upon variation of the total action in Eq. ([9]) with respect 
to X and <f>o, and setting this equal to the variation of 
the dissipation function in Eq. (TTUj) with respect to X 
and 0o j respectively, we obtain equations of motion for 
the collective coordinates 



X 2 J X' 
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(12) 



We now add thermal fluctuations that contribute as 
stochastic terms to Eqs. (fTT|) and (fl2|) . so that we obtain 
Langevin equations 



4>o{z,t) + a 



X(z,t) 



X(z,t) L z \6V eS [X,<t> ] 
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(14) 



Here we take functional derivatives of an effective po- 
tential that is a functional of X(z) and <f>o(z). Allowing 
for an arbitrary potential Vdi s (X, <f>o) due to disorder and 
inhomogeneities, this effective potential is given by 



V eS [X, 
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gB z X 



h A 



cos(20 o ) 
+ V dis (X,M , (15) 



and also includes contributions from the micromagnetic 
energy functional, the external field, and spin-transfer 



torques. Note that in the absence of current, field and 
disorder, the potential in Eq. (fT5|) is exactly the total 
micromagnetic energy E MM [^ldw} = V c g[X,(j) ]. 

The noise in Eqs. (JTHJ) and obeys 
(r)i(z,t)i~ij(z' ,t')) = aSij6(t — t')S(z — z') and 
(r]i(z,t)) = 0, where {i,j} & {X, (f> }. The strength a 
can be determined from the Fokker-Planck equation, 
which for the Langevin equations (|13[) and (|14p is given 
by22 



{1 + 0?) ?p\x,M 
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By demanding that the equilibrium Boltzmann distribu- 
tion function that follows from Eq. l[8|). given by 



(17) 



is a time-independent solution of the above Fokker- 
Planck equation, we find the strength of the thermal fluc- 
tuations as 



a = 2ak B TL z /hN 



(18) 



We see that the noise obeys the fluctuation-dissipation 
theorem with an effective temperature TL Z /N . The tem- 
perature is therefore effectively reduced by the magnetic- 
moment density in the domain wall line. 



B. Results 

Two of us analyzed the model in Eqs. (TT3]) and (TT4]) 
in the presence of extrinsic pinningi^i In this section, we 
focus instead on the clean situation, in which Vdi S = 
everywhere. In the T = case, the domain wall line will 
stay straight because the force on each point is exactly 
the same. Specializing to [3 — 0, we find that there is a 
critical current u s C rit = XK±/h. Below this critical cur- 
rent, the domain wall will not be able to acquire a finite 
drift velocity. This phenomenon is usually called intrin- 
sic pinning,— and it does not occur for (3 0. Above the 
critical current, the domain wall will acquire an average 
drift velocity, given by 



(x) = - 



1 



(1 + a 2 ) 



(19) 



In the presence of thermal fluctuations we can no longer 
assume that the domain wall remains straight, and we 
need to go through a rather more elaborate procedure to 
find the average drift velocity. More specifically, if there 
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are thermal fluctuations, we can differentiate between the 
flow regime above u S)Cr it, for which Eq. (|19[) still approx- 
imately applies, and the thermal regime below K s . cr it, in 
which the speed is finite, but goes with a different power 
law. 

To find the behavior in the thermally assisted regime, 
we start by rewriting the Langevin equations to find just 
one equation for (j>o(z,t). We then specialize to the case 
without an external magnetic field, with (3 = 0, and take 
1 + a 2 ~ 1. We assume that X" ~ because there is no 
potential that couples to X and find 



= 2a J ~^ - 



■ roj, \ aVs 
sm(20 o ) ^~ 



Vx - otr)<f> 



(20) 



which describes the motion of a string in a tilted wash- 
board potential, a problem that was investigated before 
in a different context by Buttiker and Landauer^ For 
Vs < w S)C rit) the string propagates by thermal activation. 
This occurs due to the formation of a nucleus, or of a 
kink-antikink pair in the string. That is, part of the 
string is moved over the potential barrier due to ther- 
mal activation. If this nucleus is large enough, the kink 
and antikink will proceed to move apart from one an- 
other, and the string propagates to the next potential 
valley. Two factors are important: the energy barrier 
AE that needs to be overcome to generate a sufficiently 
large nucleus, and the propagation velocity of the kinks 
(and antikinks). 

Let us sta rt wit h the former. It is given by 
AE/Eq — \JJ/K±_ J(d4>o t N/dz) 2 dz, where 0o,tv is a 
time-independent solution of the differential equation 



(J/K ± )(d 2 ^ N /dz 2 ) = sin^jv) 



v 8 /v StCI it, and rep- 



resents a stationary configuration, corresponding to the 
motion from a local maximum and back. In the above, 
AE is given in units of Eq = a^/JK±/(2X). The result 
is shown in Fig. [TJ In the limit that we are close to the 
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FIG. 1: (Color online) The energy barrier AE/Eo necessary 
to generate a nucleus large enough for the kink and antikink 
to propagate as a function of u s /« s , C rit- The solid line was 
found by numerical calculation. The dashed line corresponds 
to Eq. pTj) . 



critical current, i.e. that (t> s ,crit — v s )/v s 
solve for AE exactly and find2& 



< 1, we can 



AE 



24 s/JKl 



1 - 



^s.crit 



5/4 



(21) 



This formula indeed fits very well to the tail of our nu- 
merical curve. Note that the energy barrier remains finite 
as v s — > 0. The difference between the above formula and 
the numerical solution is at most 25%, suggesting we may 
use Eq. (f2"Tj) to estimate the qualitative behavior of do- 
main wall motion even at lower values of w s /w s .crit- In 
the limit that v s — > we have that AE cx v s log v s . This 
limit is not shown in Fig. [1] as it applies only for v s very 
close to zero. 

The other important quantity is the velocity at which 
the kink and antikink move away from one another. This 
velocity is found numerically by solving the equation 



— 4> - sin 0o 
u 



^s.crit 



(22) 



with uq = 2XEo/h, and finding the u for which 
the solution <fio(z) which starts out at 0(0) = 
arcsin(u s /w SiCr it) will go away from that point, and re- 
turn to arcsin(t! s /w SiCr it) at z — > oo. The universal curve 
for this velocity is shown in Fig [51 
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FIG. 2: Velocity of the kink (equal to minus the velocity of 
the antikink), as a function of current. 



The probability of creating sufficiently large nuclei fol- 
lows an Arrhenius law j cx exp[— AE/k^T)]. We now 
have all the necessary ingredients to find the average ve- 
locity of the string, which is proportional to \fuj. In the 
limiting case where the current is close to the critical one, 
we have that*2£ 
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(23) 



where u is a function of w s /w s ,crit as in Fig. O In Fig. [31 
we have plotted this velocity for different temperatures. 
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FIG. 3: (Color online) Velocity of the domain wall line as a 
function of current, for different values of T. 



C. Experimental Status 

In many experiments, the nanostrip is sufficiently nar- 
row that we can neglect deformations of the domain wall 
line and approximate it as being rigid, an approximation 
we treat in the next section. However, Yamanouchi et 
al^ have observed in their experiments with magnetic 
semiconductors that the domain wall looks wedge-shaped 
in the current-induced case, suggesting that deformations 
play a role in wide enough nanostrips. 

Yamanouchi et al. also found that the velocity of the 
domain wall obeys a scaling law. More specifically, they 
fitted their data with a creep-like scaling law log AT oc 
with an exponent ji ~ 0.33. Recent experiments with 
ferromagnetic metals^ have found the exponent /x ~ 1/4 
which would imply 21 that the dissipative spin-transfer 
torque dominates in the creep regime in this case. 



A. Model 

As mentioned before, rigid domain walls obey X' = 
4>' = 0, i.e., they are rigid in the z direction. We ex- 
pect this approximation to hold in the limit when L z is 
comparable to A. We integrate Eqs. (fT3j) and fT4"]) over z 
(which is quite trivial since there are no z dependences 
anymore) to obtain the Langevin equations for a rigid do- 



main wall (We use J dz6V[X(z)]/SX(z) 
and equivalently for 
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(25) 



We integrate the total potential in Eq. (fT5")) to find (note 
that X' = 4>' = 0) 



^rigid 



-hN 



k± , . ,^gB z x 

— cos(20 o ) + — - 



+ Vd is (X,0 o ). (26) 



The stochastic correlations are found from 



(f)i(t)Vj(t')) = j?(l dz I dz'rjiiz^jiz'^)) 



hN 



SijSit-t'). (27) 



In our model, rigid domain-walls obey the fluctuation- 
dissipation theorem with effective temperature T/N, i.e., 
the temperature is effectively reduced by the number of 
magnetic moments in the domain wall. 



B. Results 



Clean system, intrinsic pinning 



III. RIGID DOMAIN WALLS 

In this section, we simplify the model of domain-wall 
lines by assuming that the domain-wall coordinates are 
constant along the z direction, the domain wall is then 
rigid. This simplification allows us to obtain various re- 
sults analytically. 



We first focus on the case that the extrinsic pinning is 
zero. Substitution of X from Eq. ([24]) into Eq. ([25]) then 
provides us with an equation that is independent of X. 
Using the equilibrium solution of Eq. (|16p . we find the 
average velocity of the chirality (<fio)- With this result 
and Eq. ([24]) . we compute average drift velocities 



(X) 



k) + 13- 



gBz 
h 



(28) 
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where the average chirality velocity is given by-=i (we omit 
a factor 1 + a 2 ~ 1) 



-2ir(e H - lt -l)ak B T/hN 



(29) 



In this expression, the dimensionless effective potential is 
given by 



eri/cff 00 #1 ,„ , s 



(30) 



and the dimensionless effective field is defined as i? c ff = 
47rL z [(a — f3)v s /\ — gB z /h]/a. The expressions in 
Eqs. ([28]) and ([29]) generalize the results by Duine ei 
to include external magnetic fields and j3 ^ 0. 

In the field-driven case, we set v s = to find the be- 
havior in Fig. [U In the calculations, we use a fixed value 
for the damping parameter a — 0.02 and several values 
for the temperature. At zero temperature, the drift ve- 
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FIG. 4: (Color online) Field-driven domain-wall motion for 
several temperatures. We set v s = and take a = 0.02. 

locity depends on the external magnetic field linearly as 
(X) = X(gB z /h)/a up to a critical value H CI = aK±/h. 
At this point, the domain wall starts precessing, causing 
the Walker-breakdown. This behavior was originally pre- 
dicted by Schryer and Walker^ and was subsequently 
observed, e.g. by Beach et al~ From Fig. [4] we see 
that for nonzero temperatures, the Walker breakdown 
smoothens out. For T — > oo, it fully disappears, and we 
find that the domain-wall velocity is linear with the field 
for all fields. We note that temperature only has an ef- 
fect on the drift velocity for small fields. For very large 
external fields, the drift velocity is for all temperatures 
linear with the field and obeys (X) = a\(gB z /h). 

For the purely current-driven case, where we set B z = 
0, the relative values of a and (3 determine the sign of 
the contribution due to the current to the effective force. 
Again, we set a = 0.02 and choose several values for (3 
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FIG. 5: (Color online) Current-driven domain-wall motion for 
several temperatures and several values of j3. We set B z = 
and take a = 0.02. 



and for the temperature. It is indeed seen in Fig. [5] that 
the behavior of the average drift velocity for (3 < a is very 
different from (3 > a. In the limiting case (3 = 0, we see 
that there is a critical current w s . C rit = K±X/H, in agree- 
ment with Eq. pop . In the large (3 limit, we recognize 
a Walker-breakdown-like behavior, just like the behav- 
ior found in purely field-driven domain-wall motion. For 
very large currents, the drift velocity acquires a linear 
dependence on the spin current (X) — v s for all values of 
(3 and for all temperatures. 



2. Extrinsic pinning 

Domain- wall pinning is of practical interest. Control 
of domain walls is achieved by pinning the domain wall 
in a position by means of, for example, a deformation 
in the material. In nanowires or thin metal strips, dents 
in the sample act as an intended pinning potential. A 
thorough understanding of depinning times is especially 
important in systems that are relevant for technological 
applications, such as data-storage devices. 

To incorporate pinning in our theory, we add a pinning 
potential V p i n to the potential V r [ g id in Eqs. (f2"4"]) and (|25[) . 
A pinning potential that is due to a deformation in the 
material can be obtained through energy analysis. We 
discern between two types of deformations: symmetric 
and asymmetric notches. The sample we have in mind is 
a thin strip of ferromagnetic or semiconductor material, 
where the notches are little dents on the sides of the strip. 
In the case of symmetric notches, there is a dent on both 
sides of the sample. If there is only a dent on one side, we 
have an asymmetric notch. It turns out that symmetric 
pinning sites can effectively be described by a quadratic- 
well potential in the X direction^ that is independent 
of the chirality. For asymmetric notches, the pinning 
potential also has a chirality dependence^ 

As an example, we consider the symmetric notch. If 
we add a X dependent symmetric-notch contribution to 
Eqs. (j2"4")) and ([23)) . we see that the pinning contribution 
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in Eq. ([25]) vanishes. Because of the explicit X depen- 
dence in Eq. ([M]) . the system can no longer be described 
by a probability distribution that depends on the variable 
4>o only. Therefore, all terms in the Fokker-Planck equa- 
tion (|16|) need to be taken into account, and we are not 
able to find an analytic solution to the full problem. We 
can, however, compute depinning times using Kramer's 
escape-rate theory. 

The potential that Tatara et al£± proposed for a sym- 
metric notch has a kink at the sides, which makes it less 
suitable for escape-rate computations. We therefore use 
a very similar, but smooth potential 



\ : 



NV 



(31) 



where Vq/2 is the depth and 2£ is the width of the po- 
tential well. 

We compute depinning times for the pinning po- 
tential in Eq. (|3Tj) using Kramer's escape-rate theory 
which states that the depinning time is proportional to 
log(rwo) = AV/hkBT, where ujo oc Vo/K is an attempt 
frequency and AV is the height of the potential barrier 
that has to be overcome. We determine the positions of 
the potential minimum and the saddle through variation 
of the potential y r ; g id + Vp- ln . Substitution of these co- 
ordinates provides us with the potential difference AV. 
The depinning time for driving current and field is found 
as logr oc AV oc -{gB z /h + (3v s /\). Note that if (3 = 0, 
the depinning time is independent of the applied current 
within this approximation. 



3. Disorder potential 

The final potential that we consider is a disorder po- 
tential. Due to, for example, roughness on the edges of a 
sample, there is a random pinning potential Vdi s (X) that 
is felt by the domain wall. We consider a situation of 
strong pinning of the angle 4>o such that (<f>o) = (See 
Ref. [231 ] for an extensive analysis including the dynam- 
ics of this angle). This would be the case for fields below 
Walker breakdown and currents below u SjC rit- The equa- 
tion of motion for the coordinate X is then found from 
Eq. J25) to be 



X _ g v » 9B Z X dV dis (X) 
a A A h hN dX 



+ Vx ■ 



(32) 



The disorder potential Vdi s (X) that enters the equations 
of motion is characterized by certain spatial correlations 
[KubOX) - V dis (X')} 2 = A\{X-X')/\\-r. In this expres- 
sion, the line denotes an average over the disorder, and 
A > is a measure for the strength of the disorder po- 
tential. The exponent 7 characterizes the nature of the 
correlations. In general, there are two ways of obtaining 
a disorder potential: by applying a random field, or by 
randomizing locally the strength of certain coupling con- 
stants in the system, for example the anisotropy or the 



spin stiffness. The former is called random field disorder 
and gives rise to correlations with 7=1, whereas the lat- 
ter is called random bond disorder with exponent 7 = 0. 
For both limits, Le Doussal and Vinokur— have obtained 
expressions for depinning times for a zero-dimensional 
object (in our one-dimensional model, the domain wall 
itself is described as a point at position {X, </>o} with di- 
mension zero). 

For random field disorder (7 = 1), Le Doussal and 
Vinokur find that the drift velocity is zero up to some 
critical driving current and/or field F c cx A/T, and is 
then linear with the driving force 



(33) 



(X) _ (3v s gB z 

A A h 

For random bond disorder (7 = 0) they find that the 
drift velocity obeys 

(Pvs/x + gBz/ny- 1 



(X) 



(34) 



where z = 2 + A/2(k-QT) 2 is the so-called dynamical ex- 
ponent and r(z) is the gamma function. Surprisingly, 
they find that for < 7 < 1 the results turn out to 
resemble the results obtained for higher-dimensional ob- 
jects, i.e. Le Doussal and Vinokur find a creep scaling 
law with a certain creep exponent 



/3v s /X + gB z /H 



exp 



-(1-7) 



/3v s /X + gB z /h 
lF c 



(35) 



where F c oc T[A/2{k-B,T) 2 ] 1 ^ is a characteristic critical 
driving force. Note that the creep exponent (z — 1 — 
7/2)/(l — 7) is always larger than 1 for < 7 < 1. 



C. Comparison with other work 

To compute drift velocities at finite temperature, we 
have expanded the theory proposed by Duine et alM- to 
include external magnetic fields in addition to driving 
currents. Other theoretical work is done by Tatara et 
ali^ and by Martinez et aL— Here, we compare our re- 
sults with the results by Martinez et al. 

In order to write Eqs. ([24]) and ([25]) in terms of the 
coordinate X only, Martinez et al. assume (f>o to be small 
such that sin 20o — 20o- Note that this assumption only 
holds when B z <C gK± and/or v s <§; w s . cr it for [3 ^ 0. If, 
however, we do make this assumption, and differentiate 
Eqs. ([M)l and ([2"5|) with respect to time, we find 



X ■■ ■ K ± 
— = a(po + 2<po—— 
A n 



X 



(36) 



(37) 
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where we omitted the stochastic terms. We now substi- 
tute Eqs. dMJ) and (07]) in Eq. ^ and like Martinez et 
al. add a stochastic term to find 

n _i_ 2^_^ ^ x , R v s , gB z , 

{1 + a) 2KlJ = - a J + ^ + — + Cih - (38) 

The correlations of the stochastic force are, according to 
Martinez et al., given by 

(Cth(*)Cth(t')> = *S(t - t') , (Cth(t)) = . (39) 

From the fluctuation-dissipation theorem, they infer that 
a cx Nak-gT. An easy way to obtain the Fokker-Planck 
equation is to introduce a new variable v = X^± such 
that we have two Langevin equations 



(l + a 2 ) 



2K ± X 



X =v 



A 



A 



gB z 
h 



(40) 



(41) 



Note that the above equations have to be solved with the 
initial conditions X(t = 0) = v(t = 0) = v s to include the 
reactive spin-transfer torque. The Fokker-Planck equa- 
tion generated by these Langevin equations is given by 
(we omit a factor 1 + a 2 ~ 1) 



dP[X,v] _ 

at ~ 
_ 2K± d 

h dv 



d 

dX 



{vP[X,v}) 



v v s gB z 
^ a A +/3 A + — 



aK ± d 
h dv 



P[X,v] 



(42) 



This Fokker-Planck equation is satisfied by a Boltzmann 
equilibrium distribution that has the same potential en- 
ergy as Eq. (|26p for <j>Q — 0, but with an additional ki- 
netic energy Hv 2 /K±. The exact form of the stochas- 
tic strength is a — ak^T/NKx- We conclude that this 
procedure complies with our model for a small range of 
applicability. 

Several field-driven domain-wall motion experi- 
ments have been performed in ferromagnetic metallic 
materials ! 18 ' 32 Clear Walker-breakdown behavior is ob- 
served, but the peak is not smoothed like our theory 
predictions. Estimates by Duinc et al. 22 show that 
room temperature, at which these experiments were per- 
formed, leads to an effective temperature ksT/NK± ~ 
10 -3 for ferromagnetic metals. Our prediction is there- 
fore indistinguishable from the zero-temperature curve in 
Fig. 2J The reason for this low effective temperature is 
the fact that the number of particles in the domain wall 
is relatively high in a ferromagnetic metal. In magnetic 
semiconductors, however, not all magnetic moments par- 
ticipate in the magnetization, reducing the number of 
magnetic moments in a domain wall by up to a factor 
^100, thereby greatly increasing the effective tempera- 
ture. From Fig. 0] we see that an effective temperature 
kBT/NK± ~ 10 _1 should be distinguishable from the 



zero-temperature curve. We predict therefore that the in- 
fluence of thermal effects on Walker breakdown should be 
observable with field-driven domain walls in clean mag- 
netic semiconductors. 

Escape time studies have been performed on narrow 
domain walls by Ravelosana et al^ They indeed find 
that the logarithm of the (average) escape time decreases 
linearly with the applied current. Since the current de- 
pendence of the logarithm of the escape time is deter- 
mined by /?, we can estimate the value of (3 from their 
curve, and we find that it is of the order j3 ~ 1CP 2 , i.e. 
of the order of the damping parameter a, in agreement 
with theoretical expectations. 



IV. VORTEX DOMAIN WALLS 

In this section we consider the limit of the domain-wall 
line that corresponds to vortex domain walls and derive 
several analytical and numerical results. 

A. Model 

Vortex domain walls are described by making the 
ansatz for the z dependence of the chirality 



4> = 2arctane ±[z - z(t)1/K 



(43) 



The coordinate Z(t) plays a similar role as X(t), and 
the width of the vortex domain wall n = yJ/K±_ in 
the z direction is the equivalent of the width A in the x 
direction. The coordinates {X(t), Z(t)} now determine 
the position of the vortex at time t. 

In principle, we should consider boundary conditions 
in the z direction. We deal with this problem by as- 
suming k/L z — > 0, in which limit the ansatz reduces to 
</>o - n/2±{9[z- Z(t)} - 1/2}tt, with 9(x) the Heaviside 
step function. The drawback is that we are now neglect- 
ing all boundary effects. The ± sign is the product of 
the chirality and charge in the ansatz and determines 
whether we have a clockwise or counterclockwise rota- 
tion in the vortex domain wall, and is usually denoted 
as the Skyrmion number s e {-1,1} . We note that the 
vortex domain wall is now fully characterized by its posi- 
tion {X(t), Z(i)}, its dimensions (k x A), this Skyrmion 
number s and the number of magnetic moments in the 
vortex AnXLy/a 3 — 2nN/L z . 

We apply the simplified ansatz to the action in Eq. © 
to find an action for a vortex domain wall in terms of the 
collective coordinates X(t) and Z(t) 



Svor [X, Z] 



v<, Z , 

sn ir J ■ (44) 



The dissipation function in Eq. (|10p is also written in 
terms of the new coordinates, however, we need to take 
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into account the full ansatz in Eq. (|43[) in order to find 
the k dependence in the last term 



-Rvor [X, Z] 




2Z 2 



(45) 



Variation of the functional in Eqs. (|4"4"1) and ([43]) provides 
us with the Langevin equations for a vortex domain wall 



Z aL z X (3L z v s L z gB z 



K X 



X 



k X 



X 



2a- 



Z 



Vs v 

s^— + Vz 



(46) 



(47) 



where we have again added stochastic forces to model 
thermal effects. We see that the effective dampings are 
given by ax — ah z jn and az — 2a. We write the 
right-hand side of Eqs. ([4*6"]) and (|4T|) in terms of a total 
potential 



-kN 



0- 



,X 
A~ 



z 



gB z X 



h A 



(48) 



Note that this potential is also obtained by inserting the 
ansatz in Eq. (fTS"|) . With these identifications, we can 
again write the Langevin equations (|46|) and (|47|) in the 
more suggestive form 



Z X L z A dV vol tcx 

sir ax — = ^ 

k A kN h dX 



(49) 



TL z /Nk. The reduction by the factor kN/L z is caused 
by the fact that the number of microscopic degrees of 
freedom is proportional to this factor. Note that, be- 
cause the damping is anisotropic, it is also necessary to 
introduce anisotropy in the fluctuations. 

From the form of Eqs. (|4"9"|) and (fST))) . it is clear that 
we reach the isotropic case for ax = az, i.e. when we 
demand that L z = 2k. Note that then also ax = oz and 
that the reduction of the effective temperature is now 
proportional to N. That this is indeed the isotropic case 
is also seen from the fact that now the coordinates X/X 
and Z/k are treated on equal footing in the dissipation 
functional in Eq. (|45|) . 

If we also demand that A = re, our model describes cir- 
cular vortices, that furthermore occupy the entire width 
of the strip. This case is similar to the zero-temperature 
results of Shibata et ali^ on the current-induced vortex 
displacement in a magnetic nanodisk, and their equations 
of motion correspond to ours in the case that {3 — and 
B z = 0. 

To include extrinsic pinning in our model, we again 
add a pinning potential to the potential Vertex- As an 
example, we will consider a circularly symmetric pinning 
potential, quadratic in both X and Z ' , and bounded at a 
certain radius £ 



N 



V 



ex 2 



e)- (53) 



B. Results 



L z k <9V\ 



sir— + az — 
X K 



kN h dZ 



vortex V / cn \ 

+ Vz ■ (50) 



Note that the prefactor L z /Nk is proportional to the 
inverse of the number of magnetic moments in the vor- 
tex domain wall. Using the Fokker-Planck method out- 
lined in section II, we find that the probability distribu- 
tion function P does not satisfy Boltzmann equilibrium 
P tx exp[— Vvortex/SfcB? 1 ] if we assume that the fluctua- 
tions in the X and Z direction have the same strength. 
However, we write down a more general Fokker-Planck 
equation than the one in Eq. (flUj) . in terms of stochas- 
tic correlations (vY (t)vj (t')) = &ijd(t — t'). Again, we 
demand the Boltzmann equilibrium in Eq (|17[) to be a 
solution to the modified Fokker-Planck equation, which 
yields complicated equations that can be solved. Up to 
second order in the small parameter a, we find that the 
stochastic correlations must obey 



where the indices denote i S {X, Z} and 



o x = 2a x 



kBTL z 
HNk 



az = 2a ; 



HNk 



(51) 



(52) 



In these relations, we recognize the fluctuation- 
dissipation theorem with effective temperature T e g = 



1. Zero temperature without extrinsic pinning 

When there are no external fields or pinning poten- 
tials, we can read off the T = behavior directly from 
Eqs. (|4"6"| and (|4T)l . For instance, it is clear that under the 
influence of a current to the right, the domain wall will 
move to the right also, and its speed in the x-direction 
will be directly proportional to the current, as can be 
seen in 



X 



I + ^a(3 

2L, 



(54) 



Note that if (3 = a, the velocity of the domain wall will 
be exactly equal to the velocity of the current. For the 
transverse velocity Z, we get 



Z 



-sL z (P - a)v s 
A(tt 



77 K / 



(55) 



The direction of motion in the z-direction depends on 
the skyrmion number s, and also on the sign of /3 — a, 
and the magnitude is proportional to v s . Theoretically, 
if = a, the vortex core would move in a straight line, 
i.e. the center of the vortex would not get a transverse 
displacement. Note that there is no intrinsic pinning in 
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the case of a vortex domain wall, not even if j3 = 0. For 
/3 = a = 0, (X) = -v s , showing that the shape of the 
vortex domain wall has an influence on the motion as 
well. 



2. Zero temperature with extrinsic pinning 

We will now investigate what happens if there is a pin- 
ning potential of the form in Eq. (|53)l . If the poten- 
tial is not bounded (i.e. if the step function is absent), 
T = and if e xt = 0, the solutions can be found an- 
alytically, and they describe a circular motion ending 
at a fixed point. As we can read from the formulas, 
these fixed points will be at X ~ (3v s h^ 2 / (\ 2 Vq) and 
Z = sirvgh^ 2 / (XL z Vq) . Note that only the former de- 
pends on f3. This agrees with our physical intuition that 
(3 tilts the potential landscape in the A-direction, but 
that it has no effect in the Z-direction. 

When the potential is bounded there is a critical cur- 
rent for depinning the domain wall. We cannot find this 
current precisely using analytic calculations, but we can 
however make a rough estimation if we take the cur- 
rent for which the equilibrium position falls outside the 
boundary as an indication. This will naturally overesti- 
mate the critical current since the domain wall precesses 
after the current is switched on, but it is a good approx- 
imation for the upper boundary. We find 



^s.crit < 



\ 2 V L Z 



L 2 z /3 2 + A 2 tt 2 



(56) 



Numerical calculations can give us more precise results. 
For instance, for a symmetric domain wall (i.e. k = A) 
with L z = 100A, £ — A, and (3 — 0, we find w s , C rit — 
21.78AVb/fi., while the analytical upper bound found with 
the formula above is w s . cr it < 31.831AVo/^- Numerical 
calculations also give us the escape time at u SlC ritj which 
is approximately 0.074ft/Vb (for /3 = 0). Eq. ([51))) suggests 
that the critical current is proportional to 1//3. Fig. [6] 
shows both Eq. ([56]) and the numerical results for the 
aforementioned values of n, L z and £. 
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FIG. 6: (Color online) Critical current as a function of the 
dissipative spin-transfer torque parameter j3. 
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FIG. 7: A possible trajectory of a vortex domain wall in 
the presence of the pinning potential given by Eq. (|53[l with 
k B TL z = NVk and k = A = L z /100. The black dot denotes 
the coordinates at which the vortex leaves the potential well, 
in this example at time t ~ 1.4ft/Vb- 



3. Escape rate at finite temperature 

We now include thermal fluctuations. To give some 
insight into the motion of the domain wall in this case, 
we have plotted a possible solution to the Langevin 
Eqs. (gH|) and ([50]) in Fig. [7] The current here is just 
under the critical one, k^TL z = NVk, and we are look- 
ing at a symmetric vortex with L z = 100A. 

We next determine the escape time (i.e. the time it 
takes for the domain wall to move outside the boundary 
of the pinning potential) as a function of the current, and 
its dependence on temperature. As the average escape 
time is rather hard to determine and as the escape times 



are distributed approximately exponentially, we have, for 
practical purposes, chosen the median escape time as an 
indicator. Median/ In 2 is then expected to be a good 
measure for the average escape time. The results, with 
the logarithm of the median escape time plotted against 
the current, are shown in Fig. [8] 

Note that there is a clear distinction between the ther- 
mal regime and the slide regime. The behavior in the 
thermal regime can be fitted very well by an equation of 
the form exp[(a — bv s / v s ,c?it)NV k / k^T L z ], with a and b 
two numerical factors. Fitting shows a to be of the order 
of 35, and b around 1.8. Hence we find that the logarithm 
of the escape time is proportional to current. 
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FIG. 8: Plot of the median escape time as a function of cur- 
rent, for three different values of the temperature. 



C. Comparison with other work 

Vortex domain walls are of great interest to experi- 
mentalists, because they are large enough to be visible 
with, for instance, scanning electron microscopy or mag- 
netic force microscopy, and because their dynamics take 
place on observable time scales. Because many experi- 
ments are done with vortex domain walls, many theoret- 
ical models have been developed to describe them. He, Li 
and Zhang 1 - developed a 2D model for the vortex domain 
wall, while Shibata et al^ described the current-induced 
motion of a circular magnetic vortex. The biggest differ- 
ence between these two models and ours is that in their 
case, the vortex (domain wall) is not straight, but the 
equations of motion they find are remarkably similar to 
ours in the isotropic, symmetric case. The former group 
found the motion of the domain wall towards the edge of 
the nanostrip that our model also yields, and the latter 
paper predicted precession of the vortex domain wall in a 
potential. The model of Kriiger et al^ also predicts pre- 
cession, of an elliptical shape, and our zero-temperature 
results are in agreement with theirs. 

Current-driven vortex domain wall motion has indeed 
been observed by several experimental groupsj 20 ' 34 i 35 
Furthermore, Heyne et air* found, as expected, that the 
sign of the displacement in the Z-direction is determined 
by the skyrmion number and by the sign of (a — [3). In 
the presence of pinning, precession of a magnetic vortex 
has indeed been found experimentally as well. 19 

No detailed experimental results on current-driven 
magnetic vortices at nonzero temperature have been re- 
ported. 

V. CONCLUSIONS 

We have presented a model for the driven motion of a 
domain-wall line at nonzero temperature, and analyzed 



this model within several approximations. 

First, we considered a general domain- wall line, which 
is described as a string in a tilted washboard potential. 
We computed the average drift velocity as a function of 
an applied current for currents lower than the critical 
current. 

In the limit of rigid domain walls, we were able to 
find analytical expressions for drift velocities in the pres- 
ence of thermal fluctuations. For the held-driven case, 
the well-known Walker-breakdown behavior smoothes for 
increasing temperature. In the current-driven case, the 
drift velocity depends heavily on the ratio of the dissi- 
pative spin-transfer torque parameter [3 and the Gilbert 
damping a. Here, the curves also smoothen with increas- 
ing temperature. As a result, we found no critical current 
for (3 — at nonzero temperature. We also considered 
extrinsic pinning due to, for example, a notch in the sam- 
ple. We found that the escape time is proportional to the 
exponent of the applied magnetic field and the dissipa- 
tive spin-transfer torque. Comparison with experiment 
enabled us to estimate (3 ~ 10~ 2 for that experiment^ 3 -. 
Finally, we discussed the effect of a disorder potential on 
the dynamics of the rigid domain wall. 

For vortex domain walls, we computed domain- wall ve- 
locities at zero temperature. In the presence of an extrin- 
sic pinning potential, we found an analytic upper bound 
for the critical current. Numerical computation revealed 
critical currents just under this upper bound. At finite 
temperature, numerical simulations provided us with de- 
pinning times as a function of the applied current. We 
found two distinct regimes: one where thermal fluctua- 
tions dominate, and one where the current dominates. 
In the thermal regime, we found that the depinning time 
goes as logr cx v s /T. 

The models and results presented in this paper pro- 
vide a simple framework for describing domain walls at 
nonzero temperature. Moreover, they are easily adapted 
to situations not discussed in this paper, such as differ- 
ent geometries. We hope that our results are confirmed 
with more experimental and numerical results in the near 
future. 
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